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Monte Carlo models of DIS e7 scattering must describe hard emissions as well as 
the soft and coUinear limits. Comparison with the observed experimental hadronic 
energy flow has shown that various models underestimate the high-py contribution, 
in particular at low x. We have attempted to tune the Herwig, Pythia and 
Ariadne models to improve the agreement with the data, and to understand the 
physical implications of the changes required. Unless the physics of these processes 
is understood it will not be possible to unfold the photon structure functions 
from the e+e" data without large model dependent systematic errors. 

1 The Problem 

The measurement of F2 in deep inelastic 67 scattering, where only one of the 
electrons is "tagged" in the detector and the other one escapes unseen, involves 
the determination of the 7*7 invariant mass W from the hadronic final state. 
Because of the non-uniform detec- 
tion efficiency and incomplete an- 
gular coverage the correlation be- 
tween Wvis and W critically de- 
pends on the modelling of the 
hadronic final state. R has been 
shown EJ that there exist serious 
discrepancies in the description of 
this hadronic final state. Fig. |l| 
shows the transverse energy out of 
the plane, defined by the tag and 
the beam. For a; vis > 0.1 all of 
the generators are adequate, but 
for a;vis < 0.1 they are mutually 
inconsistent, and in disagreement 
with the data. At high i?t,out the data show a clear excess over Herwig B 
and PythiaH, while the pointlike F2GENQ sample exceeds the data. Similar 
discrepancies are observed in the hadronic energy flow per event EJ where both 
Herwig and Pythia overestimate the energy in the forward region {\rj\ > 2.5) 
and underestimate the energy in the central region of the detector. 
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Figure 1: Transverse energy out of the tag plane. 



2 The Tools 

To study the contributions of 
the various partons, the Pythia/ 
Ariadne eI energy flow in the lab 
frame as a function of pseudora- 
pidity is plotted in Figure || for 
the quark that couples to the off- 
shell probe photon 7*, denoted 
the probe quark, and for the 
quark that couples to the quasi- 
real target photon 7, denoted the 
target quark. The total energy 
flow of all partons after giuon ra- 
diation is also shown. The direc- 
tion of the tagged electron is al- 
ways at negative rj. It is apparent 
that the hump at negative rj stems mostly from the probe quark which is scat- 
tered in the hemisphere of the tag, while the hump at positive rj originates 
mostly from the target quark in the opposite hemisphere of the struck photon. 

From comparisons of the hadronic energy flow of the data with the various 
models, it became apparent that the energy flow of the probe quark needs to 
be shifted to lower 77, corresponding to an increased transverse energy. This 
can be achieved in several ways: 

Anomalous events carry more transverse momentum than hadronic events. 
Figure ^) shows the partonic energy flow after gluon radiation. Increasing the 
fraction of anomalous to hadronic or VMD type events would have the desired 
effect, but the parion density functions (PDF) used (in this case SaSlDcl in 
Pythia and GRVD in Herwig) does not readily allow changing this ratio. 

Another way to increase the transverse energy is to allow for more gluon 
radiation. This can be achieved by augmenting the inverse transverse size of 
the remnant, /i, in the Ariadne colour dipole model, shown in Figure ^d). fj, 
is set proportional to the intrinsic kx of the struck quark on an event-by-event 
basis. For VMD events, fcr is gaussian with a width of 0.5 GeV. For anomalous 
events kx follows a power law. But even a generous increase of the /i parameter 
(/i = 10) has a relatively small effect on the partonic energy flow. 

Increasing the intrinsic transverse momentum kx of the struck photon is 
another way of directly influencing the angular distribution of the hadronic 
final state. Figure ||c) shows the energy fiow for Pythia events with default 
settings and enhanced kr- This appears to be the most promising method. 




Figure 2: Energy flows of probe and target quarks. 
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Figure 3: Different methods of increasing the energy flow in the central detector region 
using the Pythia/ Ariadne Monte Carlo generators. 



3 The Fix 

In Pythia the PDF determines wliether an 
event is generated as a VMD or an anoma- 
lous event. The intrinsic kx of the quasi-real 
photon can be controlled with parametersB. 10 
Just increasing the width of the gaussian 
distribution does not produce events that 
populate the region of high St, out low x 1 , 
observed in the data (Fig. Rl). A similar de- „ ^ 

LJ Figure 4: Transverse momentum 01 

ficiency had been observed m^the resolved t^e target quark in Pythia/Ariadne. 
photoproduction data at Zeuslj, which lead 

to the introduction of a power-like /cx-distribution of the form dk^/{k'^ + ^q), 
improving the distributions of the photon remnant, fco is a constant, for which 
0.66 GeV was used. The Pythia parameters only allow adjusting the kr 
for VMD type events, figure ^, but not for anomalous evens. To change the 
intrinsic kr of anomalous events a gaussian smearing is added in quadratureu. 

Figure ^ show the i?t,out and figure |^ the hadronic energy flows on detector 
level of Pythia with default parameter settings and with the dk'^/{k'^ + fcp) 
distribution for VMD plus a gaussian smearing of the anomalous events, com- 
pared to the OPAL data taken in 1993-1995 at y^s^e = 91 GeV. In addition 
the Ariadne distributions with enhanced gluon radiation are shown. While 
the -Et.out spectrum has been improved, it still falls short of the data in the 
tail of the distribution. The hadronic energy flow generated by the enhanced 
Pythia recreates the peak on the remnant side (positive rj) seen in the data 
at low Xvis, at the expense of a somewhat worse fit on the tag side. 

The 2-jet rates, listed in Table |l|, found with the cone algorithm, requiring 
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Figure 5: PythiA: Transverse energy out 
of the tag plane. 
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Figure 6: HerwiG: Transverse energy 
out of the tag plane. 
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Table 1: Jet rates for Pythia 



Table 2: Jet rates for IfERWiG 



a minimum Ex.jct > 3 GeV in the pseudorapidity range of |77jct| < 2 are almost 
doubled over the default version of Pythia, but are still substantially lower 
than the datallj. 

Herwig separates events dynamically into hadronic and anomalous type. 
A similar dk^/ih^ + /cq) distribution of the intrinsic transverse momentum of 
the photon can be added by hand. The results of this are shown in figures || 
and 1^. Both the i?t,out and the energy flows are greatly improved with the 
inclusion of the power-like distribution, with the exception of the peak in 
the energy flow at low a:vis - high Q^, which still falls short of the data. 

The cone algorithm 2-jet rate, listed in Table |^, is more than doubled over 
the default version of Herwig and is in agreement with the data. 



4 Conclusion 

The power-like distribution of the intrinsic transverse momentum of the struck 
photon of the form dk'^/{k'^ + kg) greatly improves the hadronic final state 
distributions of both Pythia and Herwig. This improved description of the 
data should reduce the model-dependent systematic errors in the unfolded 
result of the photon structure function F2 . More fine-tuning is required. 
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Figure 7: Pythia: Hadronic energy flow 
as function of x and • The symbols are 
the same as in Fig. 5 



Figure 8: HerwiG: Hadronic energy flow 
as function of x and Q^. The symbols are 
the same as in Fig. 6 



It should be stressed, though, that this is just an ad hoc solution which 
does not explain the origin of the discrepancies. It appears that the photon 
displays a more pointlike behaviour at low x than predicted. F2 and the 
7*7 fragmentation are not orthogonal. Changing the latter will affect the 
measurement of the former. A thorough enquiry is urgently needed. 
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